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Abstract —We study the impact of phase noise on the downlink 
performance of a multi-user multiple-input multiple-output system, 
where the base station (BS) employs a large number of transmit 
antennas M. We consider a setup where the BS employs M osc free- 
running oscillators, and M/M osc antennas are connected to each 
oscillator. For this configuration, we analyze the impact of phase 
noise on the performance of the zero-forcing (ZF), regularized ZF, 
and matched filter (MF) precoders when M and the number of 
users K are asymptotically large, while the ratio M/K — 6 is 
fixed. We analytically show that the impact of phase noise on the 
signal-to-interference-plus-noise ratio (SINR) can be quantified as 
an effective reduction in the quality of the channel state information 
available at the BS when compared to a system without phase 
noise. As a consequence, we observe that as M osc increases, the 
■SINK performance of all considered precoders degrades. On the 
other hand, the variance of the random phase variations caused by 
the BS oscillators reduces with increasing M osc . Through Monte- 
Carlo simulations, we verify our analytical results, and compare the 
performance of the precoders for different phase noise and channel 
■hoise variances. For all considered precoders, we show that when 
3 is small, the performance of the setup where all BS antennas 
are connected to a single oscillator is superior to that of the setup 
where each BS antenna has its own oscillator. However, the opposite 
is true when j3 is large and the signal-to-noise ratio at the users is 
low. 

Index Terms - Massive MIMO, linear precoding, phase noise, 
broadcast channel, random matrix theory, multi-user MIMO. 


I. Introduction 

M ASSIVE multiple-input multiple-output (MIMO) is a 
promising technology for future wireless networks (TJ, 
0- This technology deploys antenna arrays containing hundreds 
of antennas, which can be exploited to significantly enhance the 
network throughput and energy efficiency performances 0. 

In particular, employing massive antenna arrays at the base 
station (BS) is expected to provide significant array gains and 
improved spatial precoding resolution for downlink transmission 
in multi-user (MU) MIMO systems |4j. This in turn is expected 
to increase the throughput per user equipment (UE), and enable 
the support of a large number of UEs at the same time. 

It is known that as the number of BS antennas M becomes 
asymptotically large, the channels between the BS and the 
different UEs become approximately orthogonal 0 , 0 . This 
indicates that significant spatial diversity can be achieved due 
to the spatial separation between the antennas of the different 
UEs. Thus, the MU-massive-MIMO downlink channel, which is 
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a MIMO broadcast channel (BC), is inherently robust to channel 
correlation effects m. 

In general, MIMO systems suffer from MU interference during 
downlink transmission, which is mitigated by means of channel- 
aware precoding methods implemented at the BS Nonlin¬ 
ear precoding methods such as dirty-paper coding are capacity 
achieving for the MIMO BC 0. However, these precoders are 
highly complex, thereby motivating the need for computation¬ 
ally simpler methods such as linear precoders 0. For MIMO 
systems, in the asymptotic regime, where M and the number 
of UEs, K , are asymptotically large, linear precoders have been 
shown to achieve close-to-optimal performance ED-GH 

One of the early works analyzing the downlink performance of 
MIMO in the asymptotic regime is ED, where it is shown that for 
M, K -A oo, a linear growth in the sum rate with M and K can 
be achieved with the zero-forcing (ZF) precoder. However, in j9j, 
it is shown that the linear sum rate growth cannot be achieved 
for the ZF precoder when M/K = 1. This issue is overcome by 
regularized ZF (RZF) precoding, whose regularization parameter 
a can be chosen such that it optimizes the signal-to-interference- 
plus-noise ratio (SINR) or the sum mean square error (sum 
MSE). In iflOl . Iflll . it is shown that, in the asymptotic limit, the 
RZF precoder maximizes the SINR in the case of independent 
identically distributed (i.i.d.) Gaussian channels. In d2l . it is 
shown that for M, K -A oo, the matched filter (MF) precoder 
requires significantly larger numbers of antennas than the RZF 
precoder to achieve the same performance. However, EO reveals 
that the MF precoder outperforms the ZF precoder in terms of 
energy efficiency, while in terms of spectral efficiency, the ZF 
precoder performs better. In some recent works fl4l . lTT5l . the 
SINRs achieved by linear precoders in the asymptotic regime are 
derived for the case where the channels between the BS and UEs 
are correlated. In addition, imperfect channel state information 
(CSI) at the BS is assumed, and the precoders are optimized 
such that the SINR achieved at the UEs is maximized. 

In most prior works on MU-MIMO downlink transmission Gl- 
ED, it is assumed that the hardware components of the MIMO 
transceiver are ideal. However, it is now well understood that the 
performance of these systems can be severely limited by impair¬ 
ments arising from nonideal transceiver hardware components 
ED. Furthermore, implementing linear precoding methods at the 
BS mandates the availability of reliable CSI. This is challenging 
since the coherence time of the channels between the BS and its 
associated UEs is finite, and thus the BS is required to update its 
CSI regularly. Also, hardware impairments affect the CSI quality 
drastically, and the phase noise caused by noisy local oscillators 
used in the transceivers is a major contributor to this problem 
M, G3, GDI- In general, phase noise manifests itself as a 
random, time-varying phase difference between the oscillators 
connected to the antennas at the BS and the UEs im, im. Phase 
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noise causes random rotations of the transmitted data symbols, 
thereby causing performance degradation. As illustrated in G3, 
|20l , phase noise also causes partial coherency loss, i.e., the 
true channel during the data transmission period can become 
significantly different from the CSI acquired during the training 
period. This is referred to as the channel-aging phenomenon (ZD. 

It is therefore expected that phase noise at the BS and the UEs 
will present a serious challenge towards realizing the unprece¬ 
dented advantages promised by massive MIMO (D- The effect 
of phase noise on the uplink performance of a massive MIMO 
system has been analyzed in lfl6l , II171 , |20l . One of the early 
works on the downlink performance of MIMO systems in the 
presence of phase noise is (22l . where the error-vector magnitude 
degradation is analyzed. However, the number of studies on 
the impact of phase noise on the downlink performance of a 
massive MIMO system are limited. Prior work on the downlink 
performance of massive MIMO systems confirms that their 
performance can be reliably predicted using large-scale analysis, 
where M, K —> 00 m, da. This is because, even though the 
SINR depends on the instantaneous values of the channel and 
other random effects, these effects become deterministic in the 
asymptotic regime da, da. Interestingly, large-scale analysis 
of MIMO systems is accurate even for practical values of M 
and K. Hence, we expect that a similar analysis using tools 
from random matrix theory (RMT) lfl4l . l34l . l35l will provide 
new and important insights on the impact of phase noise on the 
performance of massive MIMO downlink transmissions, and this 
constitutes one of the main motivations for this work. 

In this paper, we analyze the massive MIMO downlink per¬ 
formance of linear precoding schemes, including the ZF, RZF, 
and MF precoders in the presence of oscillator phase noise. We 
consider a single-cell massive MIMO system comprising one 
BS serving multiple single-antenna UEs. We analyze a general 
setup, shown in Fig. |T] where the BS employs M osc free-running 
oscillators, and M/M osc £ Z + BS antennas are connected to 
each oscillator. We refer to this as the general oscillator (GO) 
setup. Two interesting special cases arise from this general setup. 
In the first case, all BS antennas are connected to a single 
oscillator (referred to as the common oscillator (CO) setup). In 
the second case, each BS antenna has its own oscillator (referred 
to as the distributed oscillator (DO) setup). For the considered 
setups, we obtain the following results: 

• For the GO setup, we derive the effective SINR fl~3l . (26l at a 
given UE for the RZF precoder as M, K —> 00 , while the ratio 
M/K = f3 is fixed. Then, we derive the optimal regularization 
parameter a, which maximizes the effective SINR for the RZF 
precoder. Furthermore, we derive the effective SINRs of the ZF 
and MF precoders for the GO setup, by treating them as special 
cases of the RZF precoder. 

• We show that the impact of phase noise on the SINR of 
the precoders can be quantified as an effective reduction of 
the quality of the CSI available at the BS, when compared 
to the system without phase noise. As a consequence, we 
observe that as M osc increases, the SINR performance of all 
considered precoders degrades. Specifically, the desired signal 
power decreases as M osc increases for all considered precoders. 
Furthermore, the interference power increases with increasing 
M osc for the RZF and the ZF precoders. For the MF precoder, 
the interference power is almost independent from the effect 
of oscillator phase noise. However, the variance of the random 


phase variations caused by the BS oscillators reduces as M osc 
increases. 

« We compare the performance of the precoders for different 
phase noise and additive white Gaussian noise (AWGN) vari¬ 
ances by using Monte-Carlo (MC) simulations. We show that 
the SINRs derived for the precoders are accurate for relevant 
and practical values of M and K. Furthermore, we illustrate 
how the relative performance of the precoders depends on 
M osc , M, /?, the phase noise and AWGN variances, and the 
CSI quality at the BS. 

• Finally, we compare the achievable rates of the CO and the DO 
setups via simulations. A general observation for all considered 
precoders is that the CO setup performs better than the DO 
setup when /3 is small. However, the opposite is true when /3 
is large and the signal-to-noise ratio (SNR) at the UE is low. 
The remainder of the paper is organized as follows. In Section 
El we introduce the massive MIMO system model with phase 
noise and AWGN, and review the time-division duplexed (TDD) 
transmission mode and the considered linear precoders. We 
present a large-scale analysis of the received signal, the effective 
SINR, and other analytical results in Sections EH and |lv] In 
Section [V] we discuss our analytical and simulation results. We 
summarize our key findings in Section [VI] Some useful results 
from the literature are presented in Appendix A, and the proofs 
for the main analytical results are provided in Appendices B and 
C. 

Notation: Vectors and matrices are represented by boldface 
lower-case and bold-face upper-case letters, respectively. The 
complex Gaussian distribution and the real Gaussian distribution 
with mean p and variance a 2 are denoted as CAT(p,cr 2 ), and 
A ,r (//, a 2 ), respectively. The Hermitian, conjugation, expecta¬ 
tion with respect to q i, and trace operators are denoted as 
{•} H , {•}* , {•} , and tr {•}, respectively. T?{-}, 3{-}, | • | and 

Z- are the real part, imaginary part, magnitude, and the angle 
of a complex number, respectively. The L2 norm of a vector is 
denoted by ||-||. An n-dimensional complex vector is denoted by 
C raxl , while C nxm denotes the generalization to an (n x m)- 
dimensional complex matrix. The integer space is denoted by Z. 
The diag {...} operator generates a diagonal matrix from a given 
vector, while I m denotes an M x M identity matrix. 

II. System Model 

In this section, we introduce the considered single-cell massive 
MIMO system with i.i.d. flat-fading channels, oscillator phase 
noise at the BS and the UEs, TDD operation, channel estimation 
at the BS, and linear precoding. 

A. Channel and Phase Noise Models 

We consider a single-cell system consisting of a BS that serves 
K UEs. The BS is equipped with M antennas and each UE 
is equipped with a single antenna. The channel between the 
mth antenna, m £ {1,...,M}, at the BS and the fcth UE, 
k £ {1,..., K}, is assumed to be frequency-flat Rayleigh block¬ 
fading, and its gain is denoted as ~ CAf(0,cr 2 h ), where 

(j 2 h = 1. Furthermore, h^ is the ( k 1 ?n)th entry of H £ C KxM , 
which represents the channels between the BS and the UEs. 
The coherence time of the block-fading channel is denoted by 
T c . Without loss of generality, it is assumed that the large-scale 
fading component of the channel is unity. 
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Fig. 1: The general oscillator (GO) setup, where the BS has M osc free-running 
oscillators, and M/M oac £ Z+ BS antennas are connected to each oscillator. 


We consider the GO setup where the BS employs M osc 
free-running oscillators, and M/M osc £ Z + BS antennas are 
connected to each oscillator. Considering a discrete-time Wiener 

phase noise model US, ED, in the jth symbol interval, the phase 

(k) 

noise sample at the fcth UE is denoted by ip - , and that of the 
/th oscillator at the BS is denoted by o'- 1 1251 . where 

= <4-r + A/, AJ ~ Afc(0, o-l), (1) 

= #1 + A?, a; ~ AT(0, al). (2) 

Here, k £ {1 l £ {1,..., M osc }, and cr£ and cr| 

denote the phase noise increment variances at the UE and the 

BS, respectively. Since the channel is constant within T c , and 
given that the mth BS antenna is connected to the fth oscillator 
at the BS, 9^ + <j) k> is the phase noise sample that 

impairs the link between the fcth UE and the mth BS antenna. 


B. TDD and Channel Estimation 

The TDD mode of operation is assumed, where the UEs first 
transmit orthogonal pilots to the BS in order to facilitate channel 
estimation. Upon reception of the pilots from the UEs, the BS 
forms an estimate of the channel between its antennas and the 
UEs. Exploiting channel reciprocity 12, the channel estimate is 
then used by the BS to transmit data in the downlink to the UEs. 

During channel training, the UEs transmit uplink pilot symbols 
that are orthogonal in time. Specifically, the UEs transmit their 
pilot symbols sequentially in time, meaning that when one UE is 
transmitting, the other K— 1 UEs are silent. This training scheme 
permits a simple channel estimation method at the BS, and as we 
shall see later, facilitates the large-scale analysis of the MIMO 
downlink performance [[2), lfT4l . ltT31 . B1 7| , | |2()|| . Importantly, this 
training scheme aids in capturing the channel-aging effect due 
to phase noise on the SINR achieved at the UEs. Our analysis 
of the impact of phase noise on the system performance for 
this specific training scheme qualitatively also extends to other 
training schemes—see El, where a code-orthogonal training 
scheme is also considered. 

The signal received at the BS from the fcth UE at time instant 
j = 0 can be written as 


y U ,0 — ^/Pu .k ©U. t' E/. Qj.t' + w Uj0 , 


(3) 


where perfect timing and frequency synchronization are as¬ 
sumed m, HU- In ©, y u ,o = [z/u 1 o, • • •, 2/u^o ) ] T ’ where V^o 
represents the received signal in the uplink at the mth BS 
antenna. p u ^ denotes the uplink transmit power of the fcth 

UE, ©o,fc = 


diag | 


oQW „ rp OSC) rp 


1 xM/M os 


}• 


where denotes an all-one vector of length M/M osc , 

and hfc = [h ^\..., fc^ U) ] T . co,fc denotes the pilot symbol 
transmitted by the fcth UE. w Uj o = [w^q, ..., tuj^] r , where 
v,™ ~ jV( 0, cr^,) denotes the zero-mean AWGN random vari¬ 
able (RV) at the mth receive antenna. 

Based on y u ,o, a linear MMSE channel estimate is formed for 
the fcth UE at time instant j = 0. This estimate can modeled as 
a Gauss-Markov process m, m 

ho.fc = V9o,fc©o,fchfc + v / gufcW ej fe. (4) 

Here, w e ,fe £ C Mxl represents the estimation error, and its 
entries are complex Gaussian i.i.d. RVs with zero mean and unit 
variance, ho ^ £ C Mxl in (0} is the fcth row of Ho £ C KxM , 
which contains the channel estimates of all UE channels. For 
the linear MMSE, The entries of w e ^ and ho /,, are considered 
to be statistically independent of each other [24). Without loss 
of generality, we assume that go = qo,k and qi = q\k-> and set 
the channel estimate variance as er- = goof h + qi = 1- The 

r , hk 

parameter go £ [0,1J reflects the quality of the channel estimate. 
When go = 1, a perfect channel estimate is available at the BS, 
while for go = 0, the channel estimate is completely uncorrelated 
with respect to the original channel. 


C. Downlink Transmission and Linear Precoding 

Let the data transmission on the downlink from the BS to the 
UEs commence in the symbol interval j = r, where (r < T c ). 
Here, r denotes the symbol periods that have elapsed after the 
uplink transmission of the pilot symbols from the fcth UE. The 
signal received by the fcth UE at time instant r can be expressed 
as 

Vd ,t = h r,fc 0 r,fcX T + w dik ( 5 ) 

K 

~ I ^ ' \/ Pk i SO. A:i C'T.k i T tUd,/c- (6) 

k±=l 

In (O. x r £ C Mxl is the transmit signal, and in (f6|. x r is 
written as a linear combination of the data symbols tV fc-^fci £ 
{l,...,iT}, transmitted to the K UEs. The data symbols are 
assumed to be circularly symmetric, but not necessarily Gaussian 
distributed ||27) . go,k £ C Mxl is the fcth column of the downlink 
precoding matrix G 0 £ C MxK , where G 0 = [go.i, • ■ •, go,id- 
w’d.fe ~ CA/”(0 ,ct^) denotes the AWGN RV at the fcth UE. 

In this work, we consider RZF precoding in ([6|, i.e., G 0 can 
be written as m 

Go = £ (HqH 0 + Mal M ) 1 H«P *, (7) 

where P = diag{pi,... ,pk}, Pk^k £ {1,..., K} denotes the 
power allocated to the fcth UE, and the normalization parameter 
£ is set such that the precoder satisfies the power constraint 
tr (GgGo) = 1. Note that the RZF precoder in (|7} is known 
to perform better than the other linear precoders, such as the 
MF and the ZF precoders (9), ifTTI . Furthermore, the RZF 
precoder simplifies to the ZF precoder when a —> 0, i.e., 

G 0 = £Hq ^H 0 H^ P^, and to the MF precoder when 

a -A oo, i.e., G 0 = £H?P^ d- 
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III. Large-Scale Analysis of the Received Signal 
and Achievable Rates 

In this section, we use tools from RMT to analyze the received 
signal model in (0. Specifically, we present a simplification of 
the desired signal term in 0 for the GO setup when M, K —> oo, 
while M/K = /3. Notably, we will show that in the CO and DO 
setups, the multiple-input single-output (MISO) system model 
in © can be re-written as an equivalent single-input single¬ 
output (SISO) phase noise channel including the effects of phase 
noise, AWGN, and interference [25]. Furthermore, we define the 
effective SINR, and discuss the achievable rates for the GO setup. 


A. Received Signal Model 

For the RZF precoder in 0. the received signal at the fcth UE 
in 0 becomes 


Ud,r — frfc ®r,fcGroC r Wd,k 

= \/Pfehfe ©r.fcC (h^Ho + McxIm'J ho,fc c T.k 

' -v---' 

-Isig 

+ hl& T ^ (h«H 0 + Ma. Em) ' c Trk 


+Wd,k- (8) 

In ®. we have introduced the following definitions: 
H 0 _ fe = [ho,i,..., ho,fc-i, ho,fe+i,..., ho,re], P_ fc = 
diagjpi,... ,pfc_i,pfc+i,... ,Pk}, and c T _ fe 
[c ri i,...,c Tifc _i,c Tifc+ i,...,c riK ] T . Furthermore, I sig £ C 
and I? t £ C M_ ixi denote the scaling factors associated with 
the desired symbol and the interfering symbols at the /cth 
UE, respectively. The factor I s i g is simplified in the following 
proposition. 

Proposition 1: Consider an RZF precoded downlink transmis¬ 
sion from a BS having M antennas to K single-antenna UEs 
employing TDD in the presence of oscillator phase noise. Let 
a > 0, M/K = /?,/? > 1, and go 6 [0,1]. Assume that -jH h H 
has uniformly bounded spectral norm for all M. Then, the desired 
signal factor I s i g , for M,K —> oo, can be simplified to 

I s i g = v/w®TpNete^' fe ^^ fc) ) ) ( 9) 


where 


Ipn — lim — tr{A$ r }, 

M —»oo M 


f = lim 


I M (1 T m(—a)) 


M,ir->oo y m'(-a)J2 k=1 Pk 


( 10 ) 

( 11 ) 


m{—a) = --—--(12) 


t = 


m(—a) 
i(—ot) + 1 


2 otp 


(13) 


In m A$ r = diag { lJ xM/Mosc , ..., 

e p4>l Mnsc) -4>ii m(—a) in (fl2l > is the Stieltjes 

Transform of the Marchenko-Pastur Law J32] Eqs. (1.12, 2.43)], 
and m'(-a) = ^^-\ z= - a . 

Proof: Please refer to Section B of Appendix C. ■ 


Remark 1: The terms t and £ in 0 depend on a and 3, and 
captures the channel hardening effect 0 , ®, m that results 
from the averaging of the random fading channels when RZF 
precoding is used, and M, K —»• oo. The term Tpn in CEB 
captures the effects of phase noise variations at the BS between 
the training and the data transmission phases, and is given by 


, vi osc 


(14) 


1=1 


Specifically, for the 

CO 

setup, where 

A$ t 

e j(0r-0o )\ M , and the 

DO 

setup, where 

A$ r 

diag <ye 3 ^ } ,..., e J ^ M) ^ } j, (IT4t reduces to 1171 

m M—f OO I 

gj(0T- 

CO setup 


T pn —► < 

_zi 




(15) 


Inspection of Tpn in (fl4l> and (fl5l l. reveals that ZTpn reflects 
the random phase variations caused by the oscillators at the 
BS. The variance of ZTpn decreases as M osc increases, while 
its mean is zero for all values of M osc . |Tpn| represents the 
random amplitude variations in I s i g caused by transmissions 
using distributed (asynchronous) oscillators at the BS. As M osc 
increases, the mean of |Tpn| reduces from 1 to exp 
Specifically, for the CO setup, |Tpn| = 1, while for the DO setup, 
|Tpn| = exp The variance of |Tpn| is the highest for 

M osc = 2, and decreases as M osc increases. In summary, when 
2 < M osc < oo, there are random variations in Tpn- In the DO 
setufQ, where M osc = M, Tpn hardens to a deterministic value 
that depends on r and a However, this hardening effect is not 
observed in the CO setup as the phase noise caused by the BS 
oscillator does not average out. 

From the central limit theorem (as K —> oo), and since 
the symbols c T are circularly symmetric, the interference term 
I int c T _ fc in 0 is a circularly symmetric Gaussian RV for both 
the CO and the DO setups. Furthermore, this term is uncorrelated 
(and hence independent) from the signal term. In the case where 
2 < M osc < oo, I int c T _ fc is non-Gaussian, but still uncorrelated 
from the signal term. Further analysis of I int c T _ fc is relegated to 
Section HVl and Appendices B and C. Upon applying 0. (fl5l i 
in ©, we have 


Vd ,r = y/PkQoTpvi&e^ ’ fo k + I int c T _ fe + Wd,k- (16) 

For the CO and the DO setups, the MISO system model in 0 
and 0 becomes an equivalent SISO phase noise channel in (IT6t 
ll25l . However, when 2 < M osc < oo, ( 1 1 6[ i still corresponds to 
a MISO phase noise channel, since Tpn in (IT4l) depends on the 
random phase noise variations of the multiple oscillators at the 
BS. 


B. Effective SINR and Achievable Rates 

For the CO and the DO setups, we define the effective SINR 
based on the SISO phase noise channel in (Q0 as 


SINR,. = 


llsigl 2 


(17) 


l int|| 


1 It is important to note that the results for the DO setup also hold in the case 
where M os c —> oo, while the ratio M/M osc is fixed. 
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Since the phase noise in (0 drifts symbol-by-symbol, the SINR, 
which depends on r, also varies symbol-by-symbol. In order to 
analyze the achievable rate of a given UE based on the SINR in 
we model the phase noise to be a constant within a block of 
symbols lf28ll . Therefore, the SINR which is computed for a given 
r corresponds to the SINR associated with a block of symbols, 
and can be used to determine the achievable rate |29| . Note that 
this model is implicity used in |[3j, M1 6|| , 1T71 . Furthermore, this 
model is only used to evaluate the achievable rates based on the 
SINR derived, and is not required, per se, for deriving the SINR 
and the other analytical results in this paper. The achievable rate 
computed based on this model is an upper-bound for the case 
where the SINR varies symbol-by-symbol. 

Based on the effective SINR in ( ITTI i. an upper bound for the 
achievable rate of the fcth UE for the CO and the DO setups for 
a given block of symbols (i.e., given r) is l25l 

C(SINRfc) < log 2 (1 + SINRfc). (18) 

This upper bound, which corresponds to the AWGN channel 
capacity, is generally tight for low-to-medium SINR values. 
Another upper bound for the achievable rate for the CO and the 
DO setups, which is generally tight at high SINR, was derived 
by Lapidoth et al. [30), and is given by 

C(SINR fc ) < i log 2 (27rSINRfc) - i log 2 (27reT(cr 2 + <5 pn c^) ), 

(19) 


where <5 pn = 1, when M osc = 1, and 6 pn = 0, otherwise. 
In dT9b . the second term represents the differential entropy of 
the phase noise process, + S pn ((j) T — <f> o). The 

result in l fl9) holds under the assumption that the phase-noise 
process is stationary, and has a finite differential-entropy rate. 
Combing (fT8l> and (IThl i. the achievable rate for the CO and 
the DO setups can be tightly upper-bounded as C(SINRfc) < 
min{Rate in (fT8l>. Rate in (fl9H ||25) . 

There are no results available for the achievable rates for 
the MISO phase noise channel in ®, when 2 < M osc < oo. 
However, the randomness of Tpn in this case is reminiscent of 
fading channels. Assuming ergodicity for the effective channel 
in (fl6l >. the achievable rate is written as [26l Lemma 1] 


C'(SINRfc) = log 


~ log 2 



^ _|_ E|ft|Isig| 2 \ 


( 20 ) 

( 21 ) 


where IK,,, denotes the expectation operation with respect to 
the phase noise at the BS. The accuracy of this approximation 
increases with increasing M osc . Also, I s j g and I int are not 
required to be independent. This motivates the definition of an 
effective SINR for 2 < M osc < oo as 


SINR fc 



( 22 ) 


The effective SINR in (l22t reduces to ( fl7T> for the CO and the 
DO setups. Also, note that the achievable rate computed in (|2H 
does not account for the differential entropy rates of the phase 
noise processes at the BS and the UEs. 


IV. SINR Analysis and Optimal a 
In this section, we first present the analytical results for the 
SINR achievable at a given UE for the considered precoders. 
We introduce Theorem [Q which provides the effective SINR 
for the GO setup. Then, the theorem is used to obtain the 
effective SINR when the ZF and the MF precoders are used at 
the BS. Furthermore, we analytically determine the optimal a, 
which maximizes the effective SINR at a given UE for the RZF 
precoder. 


A. SINR of the RZF Precoder 

Theorem 1: Consider an RZF precoded downlink transmission 
from a BS having M antennas to K single-antenna UEs employ¬ 
ing TDD in the presence of oscillator phase noise. Let a > 0, 
/? > l,qo £ [0,1], and SINRfc denote the effective SINR at the 
fcth UE. Then, 

SINRfc — SINR rz f fc o (23) 


almost surely, and the effective SINR associated with the fcth UE 
for the GO setup is given as 


SINR 


Tzf k ~ 


Pfcf 2 goE0|r PN | 2 


n (i - %^it P ni 2 - fraa ’ )+f 


(24) 


with 


E^ITpnI 2 — E^ 

K 

*2 = Pk ' 


tr {A<f> T } 


1 

M 


1 — e 


Mo, 


e _TC V (25) 


k i =1, 

ki 


(1 + m(—a)) 2 


(26) 


where t,m(—a),£, and Tpn are as given in (fl0t-(fT3]>. Specif¬ 
ically, E^ITpnI 2 = exp t< 7 2 ^) for the DO setup, and 

E^ITpnI 2 = 1 for the CO setup. 

Proof: Refer to Appendices B and C. ■ 

Remark 2: Theorem [T| captures the effect of phase noise on 
the SINR as an additional penalty to the quality of the channel 
estimate—when phase noise is present, the quality of the channel 
estimate, qo, degrades to goE^jTpNl 2 in the GO setup (1251) . 
The quality of the effective channel estimate decreases as M osc 
increases [25|. Also, the effective quality is reduced when r or 
cr| increase. In the DO setup, the quality of the channel estimate 

diminishes by a factor of exp . However, in the CO 

setup, there is no reduction in the quality of the channel estimate 
due to noisy oscillators at the BS or the UEs as |Tpn| 2 = 1. 

Clearly, the effect of channel aging on the SINR increases as 
M osc increases—the SINR in (124) decreases as M osc increases. 
This results from the degradation of the desired signal power 
by a factor E^X^nI 2 < 1 when M osc > 2, implying that the 
desired signal power decreases as M osc increases. But, the MU 
interference power increases with M osc , as can be seen in < 124b . 
This is because the CSI quality deteriorates due to phase noise 
£3, thereby reducing the interference suppression capability of 
the RZF precoder. Specifically, to sum for the CO setup, there 
is no effect of phase noise on the SINR since the desired signal 
power and the MU interference power is the same as when only 
AWGN is present. 

In the sequel, we will use the SINR result in Theorem [T| in 
order to derive the effective SINR for the ZF and MF precoders 
(Corollaries Q] and [2] respectively). 
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B. SINR of the ZF Precoder 

Corollary 1: Consider a ZF precoded transmission from a BS 
having M antennas to K single-antenna UEs employing TDD in 
the presence of oscillator phase noise. Let a — > 0, (3 > 1, qo € 
[0,1], and assume that the minimum eigenvalue of 2 jHH h is 
bounded away from zero. Then, 

SINR fc - SINR zffc o ( 27 ) 


almost surely, and the effective SINR associated with the kth UE 
in the GO setup is given by 


SINR zffc 


_PfcgoE^TpNl 2 _ 

■ft 7 (1 — loff^l^PNl 2 ) + fr 


(28) 


where 7p\, and t 2 are as given in m, (Ei), and (Ek 
Specifically, E^ITpnI 2 = exp—for the DO setup, and 
E^IXpnI 2 = 1 for the CO setup. 

Proof: When a is arbitrarily small, m{—a) in Ell is arbi¬ 
trarily large, i.e., m(—ct) 1 l32l . and m(—a ) + 1 « m(—a). 
Applying these approximations in the SINR of the RZF precoder 
in (El) . we directly obtain (El) . ■ 

As in the case of the RZF precoder, the SINR in (El l shows 
that the effect of channel aging is more pronounced for the GO 
setup when M osc > 2. In the CO setup, the performance of the 
ZF precoder is not affected at all by the phase noise. 


C. SINR of the MF Precoder 

Corollary 2: Consider a MF precoded transmission from a BS 
having M antennas, to K single-antenna UEs, employing TDD in 
the presence of oscillator phase noise. Let M/K = (3>l,q Q £ 
[0,1], and a —> oo. Then, 

SINRfc — SINR raffc o (29) 


almost surely, and the effective SINR associated with the kth UE 
in the GO setup is given by 


SINR mffc 


Mq 0 p k E^\T PN \ 2 

(<4 + !)££=! P* 


(30) 


where 7py is given in ( [Tol l. Specifically, E^TpnI 2 = 
exp for the DO setup, and E^|Tpn| 2 = 1 for the CO 

setup. 

Proof: When a is arbitrarily large, m(—a) in El is 
arbitrarily small, i.e., m{—a) -C 1 and m(—a) + l « 1. Applying 
these approximations in (f23l) . we can write 


qoPki r n(-a) 2 E,p |T PN | 2 

m'(-a) J2 p kl + crlm'(-a)J2k=iPk 

fei =1, 

fci t £k 

goPfcm(-a) 2 AfE0 |T PN | 

TO '(-a)Ef=iPfc( 1 + cr w)’ 
gop fc m(-a) 2 ME0 |T PN | 

(^-i) zL,p>( ‘+rt>’ 

MqoPkE ^ |Tpn| 

K+iS’ 


(31) 


(32) 

(33) 

(34) 


where (ED is reduced to (El l by approximating Y2 k =tPk ~ 


PknKx = {1,... ,k - l,k + 1,... ,K}, when K oo. 

feieifi 


in m, we exploit that when a is arbitrarily large, m'(—a ) « 
2m (- a > — -2j. Furthermore, in (El l, we exploit that as a-> oo, 
am(—a) 1, which yields the final result in (El l. ■ 

Remark 3: From the SINR derived for the MF precoder in 
ffl, it can be seen that the MU interference term is almost 
independent from the effect of phase noise at the BS in the GO 
setup (i.e., M osc > 1). This is unlike the effect of phase noise 
on the interference power of the RZF precoder in d24l) or the ZF 
precoder in (El i. However, the desired signal is affected by phase 
noise in the same manner as for the RZF and the ZF precoders. 

The SINR results for the RZF, ZF, and MF precoders in (El) . 
(El l, and (l30l) . respectively, reduce to the results in m, when 
phase noise is absent and only AWGN is present, and when the 
user channels are independent of each other. Alternatively, the 
results from a can be transformed to the results in (El) . (El) , 
and (El l by changing the quality of the channel estimate from 
qo to qoE^ |Tpn| 2 , where Tpn is as defined in (ITOl i. 


D. Optimal a for the RZF Precoder in the GO Setup 

Corollary 3: Consider an RZF precoded downlink transmis¬ 
sion from a BS having M antennas to K single-antenna UEs 
employing TDD in the presence of oscillator phase noise. Let 
a > 0, M/K = f3,[3 > l,g 0 € [0,1], and let SINR*, denote 
the effective SINR achieved by the /cth UE. Then the optimal a 
(denoted as a) that maximizes SINRfc in (El l is given by 


a = 


cr 2 + 1 - E^ITpnI 2 ^ 


|2 g 2 /3 


(35) 


E^ITpni 

where Tpn is as given in (1 1 Ok 

Proof: Similar to the steps followed in M, we differentiate 
(Ell with respect to a, and set the result equal to zero, which 
yields a in (El l. ■ 

Remark 4: When perfect channel estimates are available at the 
BS, i.e., qo = 1 and |Tpn| 2 = 1, which holds in the CO setup 
or when phase noise is absent, then a = a^//3, which is similar 
to the result obtained in m. The regularization parameter d 
becomes large when the effective quality of the CSI at the BS 
is poor, which happens for low values of q o or E^IXpnI 2 (i.e., 
severe phase noise at the BS). Under these conditions, the MF 
precoder becomes optimal. For asymptotically low values of cr 2 ,, 

the optimal a in (El l becomes lim d = L , which 

r - O-2-S-0 e 0I t pnF1oP 

does not correspond to the ZF precoder, since <5^0. However, as 
f3 becomes large, or when 1 — E^TpnI 2 ? 2 —> 0, the ZF precoder 
becomes optimal. 


V. Simulation Results 

In this section, the analytical results for the linear precoders 
presented in Section [IV] are verified by comparing them against 
the results obtained from MC simulations. Even though the 
analytical results are derived for M, K —> oo, in the sequel we 
observe that these results concur with those from simulations for 
finite values of M and K. 

A. Simulation Setup 

We simulate the system model specified in ([6]) using the RZF, 
ZF, and MF precoders, and numerically evaluate the effective 
SINR in (Ell . Then, the achievable rate in the downlink for 
a given UE is computed using (ED for all values of M osc , 
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Fig. 2: C(SINR^) for the optimized RZF precoder for /3 = 5, M = 50. rx^ = 
a v = 6°, and qo = 0.9. 


unless stated otherwise. Recall that this evaluation does not 
account for the differential entropy rates of the phase noise 
processes at the BS and the UE. We will evaluate the achiev¬ 
able rates for the CO and the DO setups as C'fSINR/,) = 
min{Rate in (fill) . Rate in (IT9t T f25l when exclusively compar¬ 
ing their performances in Section IV-EI Thereby, the effects of 
both SINR and the differential entropy rates on the achievable 
rates of the CO and the DO setups are taken into account. 

The system considered consists of a single cell with a BS 
having M = 50 antennas. Setting f3 = 5, the number of UEs 
served by the BS is K = 10. The channel between a BS 
antenna and a UE is drawn from an i.i.d. complex Gaussian 
distribution, i.e., hfi’ ~ CA/"(0,1). The coherence time of the 
channel is set to T c = 100 data symbol periods, thus resulting 
in an i.i.d. Rayleigh block-fading channel. The MC simulations 
are conducted for 10000 independent channel realizations. The 
phase noise is simulated as a discrete Wiener process CD, with 
increment standard deviation = 6° ED, CD. Unless 

stated otherwise, the time elapsed between the training period of 
the /, th UE and the data transmission period is set to r = K = 10 
symbol periods. Furthermore, all UEs use the same training 
power, and for downlink transmission, equal power is allocated 
by the BS to all UEs, i.e., P = j?Ik ■ The quality of the channel 
estimate is set to qo = 0.9 M- This is reasonable given that the 
UEs can choose to transmit at power levels p u j t such that the 
desired qo is attained. 

B. Verification of Analytical Results 

Fig. |3 compares the rate achieved with RZF precoded trans¬ 
mission from the BS to the UEs for different SNR values, where 
the SNR at the fcth UE is defined as jft- The value of a used in 
this simulation is evaluated using (l35l) . We see that the achievable 
rate (l2ll for the RZF precoder based on the effective SINR in 
OH concurs with the rate achieved in simulations for all the 
values of M osc considered. Furthermore, the SINR of the RZF 
precoder decreases with increasing M osc . Therefore, in terms of 
SINR degradation due to channel aging caused by phase noise, 
the CO setup is more robust than the GO setup when M osc > 2. 



Fig. 3: C(SINRfc) for the ZF precoder for ft = 5, M = 50. a^ = 6°, 

and qo = 0.9. 



Fig. 4: ClSINRfc) for the MF precoder for ft = 5, M = 50, cr^ = cr v = 6°, 
and qo = 0.9. 


Also, the performance of the GO setup for M osc = 5 is close to 
that of the DO setup. 

In Fig. [3 we compare the rate achieved when a ZF precoder 
is used for transmission from the BS to the UEs. We observe 
that our analytical results for the achievable rate using the 
SINR in (l28l> matches with those obtained by simulations for 
all the considered values of M osc . As in the case of the RZF 
precoder, the SINR performance of the ZF precoder degrades 
with increasing M osc . 

Fig. @] compares the performance achieved using the MF 
precoder in simulations with the performance that is obtained 
using the SINR given in (l30l) . The SINR performance of the CO 
setup is better than that of the GO setup when M osc > 2. Also, 
the gap in performance between the CO setup and the DO setup 
is smaller when compared to the case when the RZF and the ZF 
precoders are used. This is because the interference power for 
the MF precoder does not depend on Tpn or qo 
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M osc 


Fig. 7: C'(SINRfc) of the different precoders in the GO setup, where 1 < M osc < 
M, M/Mosc € Z+, for M = 50,/3 = 2, go = 0.9, = 0.06°, 

r = T c = 0.25 ms (10 4 symbol periods) 1361 . and different SNR values. 


6 


4 



°0 0.02 0.04 0.06 0.08 0.1 

al frad 2 l 

Fig. 6: C(SINRfc) of the different precoders in the GO setup for M = 
50, Mosc = 5, go = 0.9, and different cr^ values. The operating scenarios 
considered are the following: Scenario (a): SNR = 0 dB, (3 = 2. Scenario 
(b): SNR = 20 dB, (3 = 2 . Scenario (c): SNR = 0 dB, (3 = 5. Scenario (d): 
SNR = 20dB, (3 = 5. 

C. Optimal a and Precoder Performance Comparison 

In Fig. [5] a, which is numerically determined, is compared 
with the analytical result given in ( |35| >. In our numerical sim¬ 
ulations, we exhaustively search for a, which maximizes the 
SINR in m, when /? = 5,M = 5 O,ct 0 = = 6°, and 

go = 0.9. Clearly, the values of a obtained from (135b agree with 
those obtained from simulations. Moreover, at high SNR, we 
observe that the optimal linear precoder is not the ZF precoder. 
Furthermore, since the effective CSI quality in the GO setup when 
M osc > 2, is lower than that in the CO setup, a is relatively larger 
for the former case. 

We compare C'(SINRfc) of the considered precoders for dif¬ 
ferent <r 2 values and operating scenarios in Fig. [6] We set 
M = 50, M osc = 5, and go = 0.9 for all scenarios. In Scenario 
(a), we set 0 = 2 and SNR = 0 dB. In this scenario, the 
performance of the MF precoder is consistently better than that 
of the ZF precoder, and approaches the performance of the 


optimized RZF precoder as a'i increases. As the SNR is increased 
to 20 dB in Scenario (b), the ZF precoder performs better than 
the MF precoder for low values of <r|. But for large <r 2 values, 
the order is reversed. 

In Scenario (c), we set the SNR to 0 dB, and increase 0 to 
5, thereby reducing the interference. In this scenario, the ZF 
precoder is seen to perform better than the MF precoder, except 
when the phase noise is severe. Finally, in Scenario (d), the 
SNR is increased to 20 dB, and the ZF precoder significantly 
outperforms the MF precoder for all <r 2 values considered. In 
summary, we conclude that the relative performance of the MF 
and ZF precoders depends on the operating scenario, which 
depends on the values of M osc , M, go, cr 2 , SNR, and 0. 

D. Example Based on LTE Specifications 

We analyze an example based on long-term evolution (LTE) 
system specifications El, El, where we account for practical 
values of T c , symbol time T s , bandwidth BW, r, center fre¬ 
quency / c , doppler spread /<j, <j 0 , and a ( „. We choose T s = 
0.032 ps, BW = 20 MHz, f c = 2 GHz, and f d = 1000 Hz 
arising from a relative velocity of 500 km/h between the BS and 
the UEs. Letting T c = 1/4/d, the coherence time is T c = 0.25 
ms. We also consider that the time elapsed between the training 
and the data transmission phase is equal to the coherence time 
of the channel, i.e., r =T C = 0.25 ms l36l pg. 99]. 

Next, we compute er 2 and cr 2 based on a Si CMOS oscillator 
technology in ll37l . Specifically, we consider an oscillator, whose 
offset level at 90 MHz is —156 dBc/Hz. This renders <r| = 
cr 2 = 10 -6 rad 2 , or 0 ^ = 0 ^ = 0.06° using |38] Eq. (4)], 
implying that high-quality oscillators are used at the BS and the 
UE. In Fig. 0 we plot the performance of the precoders versus 
M osc for t = T c = 0.25 ms (10 4 symbol periods), go = 0.9, 
0 = 2, and M = 50 for different SNR values. At SNR = 20 dB, 
we observe that the performance of the RZF and ZF precoders 
decreases by around 0.3 bits per channel use (bpcu), as M osc 
increases, and for M osc > 10 oscillators at the BS, the additional 
degradation in performance becomes negligible. One the other 
hand, the degradation in the performance of the MF precoder 
is negligible for all values of M osc considered. Furthermore, at 


I-Optimized RZF 

-ZF 

- -MF 
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SNR = 0 dB, the performance degradation for all considered 
precoders is negligible as M osc is increased from 1 to M. 

E. Rate Comparisons for CO and DO Setups 

We compare the performance of the CO setup with that of 
the DO setup by computing the achievable rate as C(SINRfc) = 
min{Rate in (fl8l). Rate in (fl9U. We set q 0 = 0.9, cr^ = — 

6°, and t = K = 25, and analyze the rate performance of the 
RZF and MF precoders for different values of f3 in Figs. [8] and 

m 

We first consider the optimized RZF precoder in Fig. [8] For 
SNR = 40 dB, the performance of the CO setup is consistently 
better than that of the DO setup as the SINR used in (fl8l> and 
(fl9l ) is large, and exclusively determines the achievable rate. In 
particular, when (3 is small, the rate in (fl8l> is a tighter upper 
bound, and the achievable rate is determined by the SINR term, 
which is relatively larger for the CO setup. As 8 increases, the 
rate in ( IT9l ) becomes a tighter upper bound, and the SINR term 
is much larger than the differential entropy term. Hence, the CO 
setup still performs better. Now, consider the low SNR scenario 
(SNR = 0 dB). Here, once again, for small 8, the rate in (IT8l) is 
a tighter upper bound than that in ( fT~9l >. and depends on the SINR 
alone. Therefore, the CO setup performs better. However, as (8 
increases, the rate in l fl9l ) becomes a tighter upper bound, and 
the differential entropy term becomes significant compared to the 
SINR term. In particular, the DO setup has a lower differential 
entropy rate as it is only impaired by the phase noise at the UE. 
Consequently, the DO setup performs better. 

A similar performance order is observed in Fig. [9] for the 
MF precoder when SNR = 0 dB. For SNR = 40 dB and low 
8 values, as before, the achievable rate depends on the SINR 
alone as the rate in (ITSl) is a tighter upper bound. However, as [8 
increases, the rate in ( fl9l ) becomes a tighter upper bound. Unlike 
in the case of the RZF precoder, now the differential entropy 
term is significant compared to the SINR term in (1 1 9| ). This is 
because the SINR of the MF precoder is significantly lower than 
that of the RZF precoder. Moreover, the difference in the SINR 
for the CO and the DO setups is not as significant as when 
RZF precoding is used, since the interference power for the MF 
precoder is the same for both the setups ( l30l ). Therefore, the DO 
setup performs better than the CO setup. 

VI. Conclusions 

In this work, we derived the effective SINR of the RZF, ZF, 
and MF precoders for the GO setup in the presence of phase 
noise. We showed that the effect of phase noise on the SINR 
can be expressed as an effective reduction in the CSI quality 
available at the BS. Importantly, the SINR performance of all 
considered precoders degrades as the number of oscillators, M osc , 
increases. This is because as M osc increases, the desired signal 
power decreases, and the interference power increases. However, 
for the MF precoder, the interference power is almost independent 
of the phase noise. Furthermore, we showed that the variance 
of the random phase variations caused by the BS oscillators 
reduces with increasing M osc . By simulations, we demonstrated 
that the SINR approximations obtained are tight, and agree with 
those obtained from simulations with remarkable accuracy for 
interesting, and practical values of M and K. 

We showed that the optimized RZF has superior SINR per¬ 
formance compared to the ZF and the MF precoders for all 



Fig. 8: C'fSINRU) = rninjRate in <18L Rate in <191 1 of the optimized RZF 
precoder for the CO and DO setups, where qo = 0.9, = 6°. and 

T = K = 25. 



Fig. 9: C(SINRj.) = minjRate in <181 . Rate in (HI of the MF precoder for 
the CO and DO setups, where qo = 0.9, cr^ = = 6°, and r = K = 25 . 

considered scenarios. Moreover, the ZF precoder performs better 
than the MF precoder when the CSI available at the BS is reliable, 
the phase noise at the BS is not severe, and (3 is large. However, 
for low SNR, severe phase noise at the BS, and small /3, the 
MF performs better. Finally, we observed that for all considered 
precoders, the CO setup has a higher achievable rate than the DO 
setup when [8 is small, while the DO setup outperforms the CO 
setup when the SNR at the UE is low and 8 is large. 

APPENDIX A 

Important Results from Literature 

Definition 1 (Stieltjes Transform J I32I Section 2.2]): Let X be 
a real-valued RV with distribution F. The the Stieltjes transform 
m(z) of F, for z £ C such that Sjz} > 0, is defined as 

m(z) = E 

Lemma 1: Stieltjes Transform of the Marchenko-Pastur Law 
[32\ Eqs. (1.12, 2.43)]: Let H £ C iCxM , with entries that 
are zero-mean i.i.d. RVs with variance 1/M. Then the em¬ 
pirical distribution of the eigenvalues of H h H converges to 
the Marchenko-Pastur law almost surely as M, K —> oo, with 


A- 


-dF{ A). 


(36) 
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M/K = (3. Furthermore, the Stieltjes transform m(z), with 
complex argument 2 : such that S{z} > 0, of the Marchenko- 
Pastur density law is defined as 


h H (U + ghh 


H\-l 


—trA 1 — m(—a) 
M v ' 


0, —trA 
’ M 


< 


c 


—trAU' 1 - — trA(U 
M M v 


qhh*)- 1 


M ’ 

x H U (A + g 0 xx H + giww 1 
t 2 (l + qih) M—f O O Q 
1 + 

w H U (A + goxx H + giww 
f 2 (l + qoh) M—fpo 


■ 5 2 XW 


Q 2 WX 


H\ 


H 


g2Xw H + g2WX H ) 


w 


1 + 

x H U (A + g 0 xx F 

— 92 * 1*2 M->o o 
1 + 1 1 

w H U (A + g 0 xx 

— g2*l*2 M-»0 0 
1 + *1 


0 


0 


H 


0. 


giww H + g 2 xw H + g 2 wx H ) 


5iww h + g 2 xw H + g 2 wx H ) 


, , 1 -zfi-p± \//3 2 z 2 — 2(/3 + l)zf3 + (1 — P) 2 ^ 

m{z) =- 2 ^- (37) 

Lemma 2 (Matrix Inversion Lemma /55] Eq. (2.2)]): For an 
invertible matrix U £ C MxM , h £ C Mxl , an£ j ^ ^ ^ where 
U + ghh H is invertible, 

h H U 1 


(38) 


1 + gh H U -1 h ’ 

since h H U _1 (U + ghh H ) = (1 + gh H U' 1 h)h H . 

Lemma 3 (Resolvent Identity fl4\ Lemma 2]): Given two in¬ 
vertible matrices U and V of size M x M, 

U' 1 - V -1 = -U _1 (U - V)V" 1 (39) 

holds. 

Lemma 4 (Trace Lemma [M\ Lemma 4]): Let x,w ~ 
CAf(0, jjl m) be mutually independent vectors of length M, 
and also independent of A £ C MxM , which has a uniformly 
bounded spectral norm for all M. Then 

x H Ax-trA —> 0, x H Aw —> 0. (40) 

M 

Lemma 5: Let H £ C KxM , M,K —y oo with M/K = f3, 
whose entries are zero-mean i.i.d. Gaussian RVs with variance 
1/M, and define A = -pH H H + oIm- Then 


Lemma 9: Let U, V £ C Mx A1 be freely independent random 
matrices li34l Page 207] with uniformly bounded spectral norm 
for all M. Further, let all the moments of the entries of U, V be 
finite. Then, 

-UrUV - -UrU-i-trV M -^° 0. (48) 

M M M 

APPENDIX B 

Extensions to Existing Lemmas in fl4ll 

First, we provide Lemma [TOl which is an extension to Lemma 
H] fl4l Lemma 7]. Lemma [lO] is then used to the derive the 
effective SINR in d22>. 

Lemma 10 (Extensions to h!4\ Lemma 7]): Consider U, A, 
N £ C AIxM with uniformly bounded spectral norms for all 
M, such that A is invertible, and N is unitary, NN H = I m- 
Furthermore, assume that N is freely independent of U, A f34l 
Page 207]. Let x,w ~ CJ\f(0, jjIm) be mutually independent 
vectors of length M, and also independent of U A. Define 
90 , 91,92 S K+ such that q 0 qi = gf, 9o + 91 = 1 and 
h = -jgtrA' 1 , f 2 = -pTrUA -1 . Then, 

x h NU (A + goxx H + 9 iww h + g 2 xw H + g 2 wx H ) 1 N H x 


~ t 2 - 


go*i*2 
1 4 -1 \ 


tr{N} 


M 


M—t oo 


(49) 


x H U (A + go xxIi + giww 11 + g 2 xw H + g 2 wx H ) 1 N H x 
*2(1 + gi*i) tr {N H } m -¥0 0 q 
lTh M 

w H U (A + goxx H -(- giww 11 + g 2 xw H + g 2 wx H ) 

—g 2 fif 2 tr {N H } m—► oo 


— m'{—a) —> 0 (41) 


Lemma 6 (Rank-1 Perturbation Lemma 1 33] Lemma 2.6]): 

Let £ > 0, U, and A £ C MxM with U being nonnegative 
Hermitian, h £ C Mxl , and q £ K. Then, 

|trA [(U + CIm + ghh 11 )’ 1 — (U + CIm) -1 ] ' ' HA| ' 


(42) 


(50) 


Hrl N H ) 


1 + *1 


M 


0 


(51) 


Lemma 7 ( l[T4\ Lemma 6]): Let U, A £ C MxM with U 
being nonnegative Hermitian, h £ C M and q £ R, then, 


Proof: Since N is freely independent of U, A, from Lemma 
[9] we have 


(43) 


—trNA -1 - — trN — trA" 
M MM 


l M—foo 


0 


Lemma 8 ( HT4\ Lemma 7]): Consider U, A £ C MxM with 
uniformly bounded spectral norms for all M with A being 
invertible. Furthermore, let x,w ~ CAf(0, be mutually 

independent vectors of length M, and also independent of U, A. 
Define g 0 ,gi,g 2 £ R + such that g 0 gi = g|, go + gi = 1 and 
*1 = -gtrA' 1 , t 2 = -gtrUA -1 . Then, 


-J-trNUA -1 - -UrN -JrtrUA 

M MM 


~l M—f00 


0 . 


(52) 


(53) 


Furthermore, 


x h NUVx - 


trN t 2 (l + gHi) m- 


(44) 


x h A _1 N ±1 x - 


H 


M 
trN H 
M 1 


1 + fi 


M—f 00 


(45) 


(46) 


(47) 


x h UA _1 N h x - ^S-t 2 


0 

M—f 00 


0 


x h NUVw - 


M 

— 92*1*2 tl'N M-y 00 


1 +ti M 
x h NUA _1 N h x - t 2 
w h A _1 N h x M -^° 0. 


(54) 

(55) 

(56) 

(57) 

(58) 

(59) 


In order to obtain ( l54t and ( l57l >. the result in ( |53] > along with 
Lemma [8] is applied. The results in (l55l ). ( l56l >. (1581) . and 
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are obtained by using (l52t . (l53k and Lemma [4] Now, define 
V = (A + <7o xxH + 9iww h + g 2 wx H + g 2 xw H ) _1 , then 

x h NUVN h x 

= x h NUA _1 N h x - x h NU(V _1 - A)A^N h x 


H 


g 2 wx 


(60) 

H 


= x h NUA _1 N h x - x h NUV (g 0 xxH + ®ww‘ 

+ g 2 xw H ) A _1 N H x (61) 

= - 9o xH NUVxx h A 4 N h x - gix H NUVww H A’ 1 N H x 
- g 2 x H NUVwx H A 4 N H x - g 2 x H NUVxw H A _1 N H x(62) 


= t2 ~ 


= t2 ~ 


90*1*2(1 + qit\) 


1 + *1 


tr{N} 


M 


+ 


qltft 2 
1 + 1 \ 


tr{N} 


M 


90*1*2 
1 + *1 


tr{N} 


M 


(63) 


(64) 


In (l60t . the resolvent identity in Lemma [3] is used. Upon 
simplifying (162} by using (|54} - (|59} . the expression in (163} is 
obtained, which is further simplified to the desired result in (f64b 
using the fact that <?o9i = q%. 

Consider the term x H UVN H x, which is written as 

x h UVN h x 

= x h UA _1 N h x - <7o xH UVxx h A _1 N h x 

- 9 ix h UVww h A _1 N h x - g 2 x H UVwx H A _1 N H x 

- g 2 x H UVxw H A 4 N H x 
= x H UA _1 N H x - go xH UVxx H A _1 N H x 

- g 2 x H UVwx H A 4 N H x 

90 * 1 * 2(1 + 91 * 1 ) . tr {N H } 


= * 2 - 


1 + *1 

*2(1 + 9 i*i) * r {N H } 


1 + *i 


M 


(65) 

( 66 ) 

(67) 

( 68 ) 


1 +*i M 

In ( l65t . Lemma [3] is applied, and ( 166} is reduced to ( 167} using 
the result in < [59t . Applying Lemma [8] along with the results in 
(155} and (l56l >. the expression in (167} is obtained. Then, using the 
fact that qo9i = 9 2 , (167} reduces to the desired result in (168} . 
Finally, consider the term w H UVN H x, which reads as 

w h UVN h x 

= w h UA 4 N h x - g 0 w H UVxx H A _1 N h x 
- giw H UVww H A _1 N h x - g 2 w H UVwx H A 4 N H x 


- g 2 w H UVxw H A 4 N H x 


(69) 


= —<j 0 w H UVxx H A _1 N h x - g 2 w H UVwx H A _1 N H x(70) 
( 9o92*i*2 92*1*2(1 + 9 o*i) \ tr{N H } 


\ 1 + t\ 1+^1 

9 2 *i * 2 * r {N H } 


M 


(71) 

(72) 


1 + fi M 

The resolvent identity in Lemma [3] is applied to obtain 
which is reduced to (f70l> using ( |59} . Furthermore, Lemmas Q] [8] 
and the result in ( 155} are used to simplify ( 1711 ) to £72}. ■ 


Appendix C 
Proof of Theorem|7] 

The evaluation of the SINR in (122} involves computing the 
following terms: (i) the normalization constant (ii) the signal 
term I s j g , and (iii) the interference power ||I int || .In the sequel, 


we derive the three terms of interest. Also, when deriving the 
signal term I s i g , we prove Proposition [7] 


A. Derivation of £ 

First, we define Xo — +al m, and denote by Xo-k the 

matrix obtained upon removing the fcth column of Xu- The RZF 
precoder Go in £7} satisfies the power constraint tr (G^Go) = 1 
implying that 

f tr |PH 0 (h«H 0 + MoLm) ' H?| = 1 (73) 

K / , -2 

eYjp^Kk + Mal ^ + KA) Kk = 1 

fe =1 

(74) 


The expression in ((74} is rewritten by applying Lemma [2] twice 
as 


1 

M 




1 

M 


^ 2 E 


MPfehJfeX 0 2 _ fc ho,fe 
(l + ’M^ l O,kX 0 - k ^o,k) 
tr{%0 2 } _ 

(i + a^trjx - 1 }) 2 


£ 


M 



ptm'(-a) ' 
(1+m(—a))^ 


(75) 

(76) 

(77) 


Since ho,fc is independent of x 0 _ fc , and Xo -k is non-negative 
Hermitian, Lemmas 0] [6] and [7] are invoked in order to yield (176} 
from (1751) . Since the entries of Ho are i.i.d. complex Gaussian 
RVs, Lemma[5]is applied to £76}, finally rendering £77} following 
straightforward algebraic manipulations. 


B. Derivation of I s i g 


Define 

diag je-^o * 


x 0 ,fe 


(hnh 


0 n fc 


where $0 


1 T pJV 0 


4*0 = 


diag | f 


J osc 

d 1 ) 




AM) 


|. Specifically, 
| for the DO setup, and 


em e^o 

$0 = e^°I M for the CO setup. Then, the signal component is 
written as 

uTi 


Isig — 


M 


'-K & r,kXo h| 


- lu* 

0 ,k 


v'Pkqot, 

M X °’ fe 

y/Pkqi£ fj 

y/Pkqoi^Yl 

' M 

y/Pkqi£ fj 


H A* 


rXo lx o+e j( ^ ) V ° 5 ) 


X 0>- A$ rXo We,fe e- 




(fe) 


x 6)fc A ^rXo x 0 ,k& 




- 1 , 


M 


X 0,k A ^rXo^e,k'e 


*p]<,V < f ) -V o fc) ) 


(78) 


(79) 


(80) 


xM/M 0 . 


In £ 79 }, A$ r = diag { e^-^ol J 

e^ M ° Bc) -<*r° sc) lJ xM/Mnsc }. For the CO setup, A$ r = 

e j^-r-<h)i M) and A$ r = diag { e + , r 1) --*o 1> ,..., e ^r M) -^ 0 M) } 

for the DO setup. 

The expression in £79} is obtained by using the channel 
estimate £4} in £78}, and £79} is further rewritten as £80} by letting 
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IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY 


W e ,fc = W ejfe e 


(fc) 

, given that w 6j fe is circularly symmetric. 
Define t\ = -jtr ^ 0 1 _ A . j- - Upon invoking Lemma |9j for the 

DO setup we have -ptr |a^> t x 0 1 _ /s | = -p-tr {A<F T } U, where 

A<F r is freely independent of Xq 1 . This is because <I> r and Xq 1 
contain entries that are statistically independent of each other, 
and Xo 1 is unitarily invariant lf34l . l32tt . [|35l Theorem 22.2.3]. 
Finally, by employing Lemma [8] i.e., (144b . the signal and noise 
part in (l80b becomes 

1 H \ j. -1 -p-tl {AtF T } fi(l + Qlti) M—yoo 

M X " ' " 

1 


i x o.k A$ TXo lx o,t - —— - h J + ^~ o ( 81 ) 


T7 X OT A ^rXo W e,fc “ 


-92^tr{A$ T }f? 


M ^ '~ u c,,v h +1 

Thus, the desired signal can be written as 

\/PkQo£jfto {A$ T } fi(l + qiti) 


0, 


(82) 


Isig — 


(go + qi)ti + 1 

2 "I 


, -,-92^tr{A$ r }f : 

+\/Pkqi£,- 


,3 ('pW-'pM) 


(83) 

(84) 


ti H- 1 

= VWk£m(-a)jitr{A$ T } ^ ((p w_ (p vc) ) 
m(—a ) + 1 

where Lemma [5] is used to reduce (l83b to (l84b . and this corre¬ 
sponds to the result in ([DJ (i.e.. Proposition [TJ. 

T„l 

For the DO setup, -ptr {Ad? T } = e __ 2 ~ |fl7l . and for the CO 
setup, -p.fr {A$ t } = eA^-Vo). When 2 < M osc < oo, 

^osc 


—tr{A$ r } = —— V 
M 1 1 A/f _ _ _ Z-* 


1 




M 




(85) 


i=i 


and with straightforward calculations, it can be shown that 

2 


E,* 




1 — e 


Mo: 


( 86 ) 


C. Derivation of ||I int || 

Define Xo - -fcP-fc^o.-fcXo 1 . *2 = 


iftr {x 0 l fe HH_ fc P_ fe H 0 ,-fcXo 1 }’ and x o ,k = Then, the 

power of the interference signal I int in (El l is evaluated as 


T h T 

A int A mt 


= ^hT© Tifc Xo 1 H». fe p - fe H 0rfe Xo 1 ©; fe h fe * 


M 2 

( 2 uT 


(87) 


= ^h^rXo^H^-fcP-fcHo^Xo 1 *;^* 

e 


+ 


Jlf2^ k ^ T (*0 X 0 ,-fc) 


H 0 H _ fc P- fc H 0 ,- fe 




■Xo 

= ^h^ T x-^H». fc P_ fc H 0rfc Xo^;h fc * 

£2 _ 

- Jp^I^rXo 1 (Xo ^ Xo-fc) Xo^fc H o!-fc P -fcH 0 -A 

■ Xo^;h fe * 


( 88 ) 


(89) 


-,H A* - A,,.* . „H 


= J^^r,k^rX Q ^>r,k ~ 

■ (g 0 x Tifc x^ fc + giw* iT)fc wJ iTifc + g 2 x Tifc w^ Tfc 
+92W* iTifc x“ fc ) X 0 A $* x r,fc- (90) 


The resolvent identity lemma (i.e.. Lemma [3]) is applied in (l88b 
resulting in the expression in (f89l >. The definitions of Xo> and 
Xo.-fe, and the channel estimate in ifTJ are employed in (f89l) to 
obtain d90b . Finally, the quadratic forms in (l90l) are simplified 
using Lemmas [8] and [TO] as 


x^A$ T x 0 A$’ x T ,t qotit 2 
-f 2 H— 


tr{A<l> r } 
M 


M—> oo 


M 2 

-l. 


X r,fe A ^^X 0 X r,fc 


M 2 


1 + t\ 

+ gl^l) M—yoo 

1 + 1 1 


0 


(91) 

(92) 


XrV^X^W^t _ -g 2 *l M—yco 

M 2 1 + 1 1 ED 

X r,fcXo A ^r x T,fc ~ t ^ 2(1 + <7lfl) M—>oo q 

M 2 l + ti W 

W u,T,fcXo A ^ > r x r,fc —q 2 tlt2 M " A T—»p o „ 

M 2 l + ti ED ‘ 


(93) 

(94) 

(95) 


In order to evaluate f 2 , we invoke Lemmas [6] and [7] to obtain 


tr 


{ p -A,-i (a„A-W)h 0 h -»} 

M 


M- 


0. (96) 


Thus, for M —>■ oo, we simplify i 2 « A tr |p_ fc H 0 _ fe Xo 2 H“_ fc | 
as 


K 1 

^ ^ Af Pfcl ^°- fc 1 ^0,/ci 

(97) 

k\^k 


K ^ 

_2 hS, fcl (98) 

fc i=i, 

ki^k 


K d-h T v -2 h* 

v~^ AT 11 0,fciA.g_ fel 11 o,fc 1 

^ Pfel (l + AhT v -i L* ) 2 

fci=i, ' M u 0,fci ^0 

ki^k 

(99) 

^ CT fiig tr {Xo 2 } 

> Pfci r a 

(l + a^trjxo 1 }) 2 

ki^k ^ J 

(100) 

K K \ 

rri (—a) 

, Pkl (1 + m(-a)) 2 

(101) 


fei=i, 
ki^k 


In order to obtain ( l99l . Lemma [2] is applied twice to ( l98l >. 
Further, Lemmas |U [6] and [7] are employed in (199b to obtain 
(11 OOl i. and Lemma 0 is applied to (1 100b yielding (1101b . Hence, 
the interference power can be reduced as 


C 2 (, qotih 
My 2 1 + ti 


tr{A$ r } 


M 


e_ 

M 


( 

tr {A^? t } 

2 f 1 ^ 2(1 + gHi) 2 

tr{A$ r } 

2 g|*?* 2 

l 

M 

(l + ii) 2 

M 

(l+ii) 2 

-2 

tr{A$ r } 

M 

2 < 72 * 1 * 2(1 + giii)\ 

(l + ti ) 2 J 


(102) 
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e 

M 


( 


^2 — 


Qotlt2 


tr{A^» r } 
M 


qotit2 


tr{A«l» T } 

M 


1 + ^1 (1 + ^ l ) 2 

where t 2 is given in (1101 b . and t\ = m(—a) for M 00 . 


(103) 


V 
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